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EARTHQUAKE PREDICTION METHOD AND SYSTEM 



Technical Field 

The present invention relates to an earthquake prediction method and 
5 system, and more particularly to an earthquake prediction method and 

system for measuring the magnetic field of each point within an observation 
area to predict earthquakes. 



Background Art 

10 The Japanese archipelago is situated on volcanic chains, and the 

occurrence of a massive earthquake disaster is therefore a major concern in 
the oceans to the east, southeast and south, not to mention throughout the 
archipelago. To protect the lives and property of citizens from earthquake 
disasters, more suitable measures must be taken not only for systems for 

15 recovery following an earthquake, but for earthquake prediction, and for this 
purpose, the establishment of a low-cost, and moreover, highly accurate and 
effective earthquake prediction technology is urgently required. 

From ancient times, many examples of precursors of impending 
earthquakes have been reported, such as the agitated behavior of catfish or 

20 the increased activity of rats. In addition, several days before the Great 

Hanshin-Awaji (Kobe) Earthquake and Fire of January 1995, local amateur 
radio operators observed radio-wave abnormalities, and a phenomenon 
known as an "earthquake cloud" was witnessed by many people. 

It is believed that the occurrence of such radio-wave abnormalities 

25 and "earthquake clouds" result from some form of influence of telluric 

currents that are in turn generated by the piezoelectric effect resulting from 



the collision of tectonic plates. 

The prediction of earthquakes by observing the amount of rotation of 
a permanent magnet that is suspended from a thread has been disclosed 
(for example, in Japanese Patent Laid-Open Publication No. H1 1-258353). 
5 However, the prediction of earthquakes with high accuracy by 

observing changes in telluric current would necessitate the establishment of 
a large number of closely-spaced observation facilities over an extensive 
observation area. As an example, earthquake prediction in the eastern 
seaboard region of Japan would require an observation system covering 

10 several hundred kilometers square. 

With the current financial difficulties faced by national and municipal 
governments, such an establishment of a large number of observation 
facilities would be a tremendous burden, and the costs for maintaining and 
operating the observation system would render the realization of such a 

15 system highly problematic. 

Disclosure of the Invention 

It is an object of the present invention to provide an earthquake 

prediction method and system for collecting magnetic field data at numerous 
20 points within an observation area at low cost, and moreover, within a short 

time period to thus enable accurate earthquake prediction. 

In the earthquake prediction method of the present invention, 

magnetic fields are measured at various points within an observation area to 

both estimate telluric current induction fields and to thus estimate telluric 
25 currents; and the state of the telluric currents within the observation area and 

the change of the telluric currents over time are then analyzed to estimate a 
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seismofocal zone, a time of occurrence of a seismic event, and seismic 
intensity. 

In addition, a magnetic field noise component is eliminated from the 
observed magnetic fields, the amount of divergence of the direction of the 
5 magnetic fields from the direction of true north at observation points from 
which this magnetic field noise component has been eliminated is found, the 
telluric current induction fields are then estimated based on the vector 
difference between the observed magnetic fields from which the magnetic 
field noise component has been eliminated and a geomagnetic vector that 
10 has been corrected to true north. 

Further, the estimated telluric current induction fields are plotted on a 
map; the telluric currents are estimated by both joining points on the map at 
which geomagnetic abnormalities are noted and applying the right-handed 
screw rule; and the region in which the estimated telluric current is 
15 concentrated is estimated to be a seismofocal zone. 

Still further, past data of telluric current induction field intensities at 
specific observation points are gathered and a telluric current induction field 
intensity change pattern that indicates changes over time is generated; and 
this pattern is then compared and collated with past telluric current induction 
20 field intensity change patterns that have been stored to estimate the time of 
occurrence and seismic intensity of a seismic event. 

The earthquake prediction system of the present invention is provided 
with: mobile units such as vehicles or ships, each unit carrying a magnetic 
force line sensor for supplying as output magnetic field data that indicate the 
25 direction and intensity of lines of magnetic force, a GPS position detector for 
receiving radio waves of a GPS satellite and supplying as output position 
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data that indicate position, and a data transmitter for transmitting these data; 
and an earthquake prediction center for collecting these data of various 
points that are transmitted by the mobile units that travel within an 
observation area and for then implementing earthquake prediction. 
5 In addition, the earthquake prediction center includes: a data receiver 

for receiving data that have been transmitted from the mobile units by way of 
a communication network and an antenna; a data storage unit for retaining 
and storing various data such as map data and data that have been received 
by this data receiver; a telluric current induction field estimation unit for 

10 estimating telluric current induction fields based on map data and data that 
have been retained and stored in this data storage unit; a telluric current 
estimation unit for estimating telluric currents based on the telluric current 
induction fields that have been estimated; a telluric current induction field 
intensity change pattern generation unit for accumulating transitions over 

15 time of the telluric current induction field intensities and then generating 

change patterns; and an earthquake prediction unit for analyzing the telluric 
currents that have been estimated and the change patterns of the telluric 
current induction field intensities to estimate a seismofocal zone, seismic 
intensity, and time of occurrence of a seismic event. 

20 Alternatively, when a mobile unit is provided with a car navigation 

system, the position data of the car navigation system may be used instead 
of data from a GPS position detector. 

In addition, magnetic force line sensors and communication 
equipment may be installed in preselected existing fixed structures within the 

25 observation area, and the communication equipment may transmit the 

magnetic field data output of the magnetic force line sensors and information 
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indicating the installation positions to the earthquake prediction center by 

way of an existing communication network. 

Further, a magnetic force line sensor and GPS position detector may 

be incorporated in, for example, a mobile telephone, and the mobile 
5 telephone may use its own communication capabilities to transmit 

observation data to the earthquake prediction center. 

Still further, acceleration sensors may be provided in existing fixed 

structures such that magnetic field data are transmitted when the 

acceleration sensors detect earthquake motion; or acceleration sensors may 
10 be provided in mobile units or in mobile telephones such that magnetic field 

data are transmitted when an acceleration sensor detects a stationary state 

of at least a fixed time interval. 

According to the present invention, by installing magnetic force line 

sensors, GPS position detectors, and data transmitters in, for example, 
15 vehicles or ships that travel within an observation area, collecting at an 

earthquake prediction center the magnetic field data at each point within the 

observation area, and analyzing the telluric current induction fields and 

telluric currents that are estimated based on the magnetic field data, 

accurate earthquake prediction can be realized at a low equipment cost 
20 without necessitating the installation of measurement equipment at a 

multiplicity of points. 

In addition, by establishing magnetic force line sensors at preselected 

existing structures within an observation area and then transmitting magnetic 

data to an earthquake prediction center by way of an existing communication 
25 network, or alternatively, by establishing magnetic force line sensors in 

devices such as a mobile telephones and then transmitting magnetic data to 
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the earthquake prediction center, the observation data of a multiplicity of 
points within an observation area can be collected at a low equipment cost 
to realize accurate earthquake prediction. 

Further, the installation of a magnetic force line sensor together with 
5 an acceleration sensor whereby observation data are automatically 
transmitted when the acceleration sensor detects earthquake motion 
enables the collection of magnetic field data at the time of the occurrence of 
foreshocks that precede a main shock and thus allows observation data to 
be obtained that are effective in earthquake prediction. 

10 

Brief Description of the Drawings 

FIG. 1 is a block diagram showing an embodiment of the present 
invention; 

FIG. 2 is a view showing a model of telluric current in the vicinity of a 
15 seismofocal zone preceding the occurrence of an earthquake; 

FIG. 3 shows the earthquake prediction operation of earthquake 
prediction center 4 that is shown in FIG. 1 ; 

FIG. 4 shows the estimation of a telluric current induction field; 
FIG. 5 shows an example of the estimation of a telluric current; 
20 FIG. 6 shows an example of telluric current induction fields and 

telluric currents that have been measured; 

FIG. 7 shows an example of telluric current induction fields and 
telluric currents that have been measured; 

FIG. 8 shows an example of telluric current induction fields and 
25 telluric currents that have been measured; and 

FIG. 9 shows an example of the pattern of change in telluric current 
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induction field intensity. 

Best Mode for Carrying Out the Invention 

Explanation next regards the present invention with reference to the 
5 accompanying figures. 

FIG. 1 is a block diagram showing an embodiment of the present 
invention, and shows an earthquake prediction system that uses: mobile 
units 1 such as a vehicle or a ship that are capable of movement on land or 
on sea; or existing fixed structures 2 that have been selected in advance 
10 within the observation area, to collect observation data at a multiplicity of 
points within the observation area and predict, for example, the seismofocal 
zone and time of occurrence of an earthquake. 

In this case, the earthquake prediction system is provided with: 
mobile unit 1 in which, for example, magnetic force line sensor 1 1 and GPS 
15 position detector 1 2 are mounted and which transmits observation data to 
earthquake prediction center 4; existing fixed structure 2 in which magnetic 
force line sensor 1 1 and communication equipment 14 have been installed; 
communication network 3 for transmitting observation data to earthquake 
prediction center 4; and earthquake prediction center 4 for predicting 
20 earthquakes based on the observation data of a multiplicity of points within 
the observation area. 

Mobile unit 1 is vehicle 1-1 or ship 1-2 that moves within the 
observation area and that carries: magnetic force line sensor 1 1 that 
supplies magnetic field data that indicate the intensity and direction of a line 
25 of magnetic force; GPS position detector 1 2 that receives radio waves of a 
GPS satellite and that supplies position data; as well as data transmitter 13 
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that transmits observation data to earthquake prediction center 4. 

When the mobile unit carries a car navigation system, the position 
data of the car navigation system can be used and GPS position detector 12 
can therefore be omitted. 
5 The observation data may be transmitted in real time, or a data 

storage apparatus may be provided for recording the observation data. It 
goes without saying that the observation data may also include, in addition 
to magnetic field data and position data, data indicating the observation time. 
Setting the observation data for automatic transmission at 
10 observation positions and times that have been set in advance allows a 
reduction of the burden imposed on human operators. 

In addition, effective observation data can be obtained if the entire 
observation area is observed by mesh coverage. Mesh coverage is a 
covering method that is frequently used in radio-wave checks for mobile 
15 telephones and involves preparing a map in advance of the region that is to 
be covered, casting lines in a mesh pattern, and then traveling in a mesh 
pattern along these lines. 

Existing fixed structures 2 that can be considered include: 
electrical/gas/water meters that are established in each residence and 
20 business; vending machines that are established along roads; poles for 
power lines and communication lines; poles for traffic signals; operation 
display apparatus that is installed at bus stops; and further, the buildings of 
mobile telephone base stations and PHS base stations. 

Communication equipment 14 transmits either by cable or by radio to 
25 earthquake prediction center 4 the magnetic field data that are supplied as 
the output of magnetic force line sensor 1 1 along with information indicating 



the position of installation. 

In this case, transmission may be carried out in real time, or 
transmission may be implemented automatically at predetermined 
observation times. 

5 The data transmission path can be easily secured when already 

existing fixed structure 2 is a pole for power lines or communication lines, a 
pole for a traffic signal, or the building of a mobile telephone base station or 
PHS base station. In addition, lower costs can be obtained if an existing 
wireless communication means such as an emergency radio system is used. 

10 When an installation such as an electrical/gas/water meter is used, 

the automatic meter reading system may be used for transmission. If the 
operation display apparatus that is installed at a bus stop is used, the vehicle 
operation management system may be used for transmission. Further, if 
vending machines are used, a system may be constructed for transmitting 

15 the sales or stock information of the vending machines along with the 
magnetic field data. 

Communication network 3 is an existing communication network such 
as a mobile communication network that includes base stations or a satellite 
communication network that transmits by way of communication satellites. 

20 Earthquake prediction center 4 includes: data receiver 41 for 

receiving observation data by way of communication network 3 and an 
antenna; data storage unit 42 for retaining and storing various data such as 
observation data and map data; telluric current induction field estimation unit 
43 for estimating telluric current induction fields based on the observation 

25 data and map data that are retained and stored in data storage unit 42; 

telluric current estimation unit 44 for estimating the telluric currents based on 
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the results of estimating the telluric current induction field; telluric current 
induction field intensity change pattern generation unit 45 for accumulating 
transitions of the telluric current induction field intensity over time and 
generating change patterns; and earthquake prediction unit 46 for analyzing 
5 the estimation results for telluric currents and the change patterns of the 
telluric current induction field intensity and then predicting the seismofocal 
zone, the seismic intensity, and the time of occurrence of an earthquake. 

FIG. 2 shows a model of telluric current in the vicinity of seismofocal 
zone of an earthquake preceding the occurrence of the earthquake. 

10 In this case, plate A and plate B are moving in mutually colliding 

directions and thus press against each other. Point C at which the local 
pressure in the boundary plane between plate A and plate B is building is the 
seismofocal zone. 

At this seismofocal zone C, massive stress is concentrated and an 

15 extremely high-stress state is in effect, and this stress gradually rises with 
movement of the plates. During this state, it is assumed that voltage is 
generated by the piezoelectric effect in seismofocal zone C, and that charge 
within the rock mass flows into seismofocal zone C. 

It is further assumed that the flow of charge within the rock mass 

20 (telluric current) follows points having good conductivity within the rock mass 
and flows from each direction like rivers. Normally, since most seismofocal 
zones are underground, most of the charge flows underground and it is 
assumed that the charge rarely flows on the earth's surface. 

It is hypothesized that immediately before the rock mass of 

25 seismofocal zone C fractures, the telluric currents rise with increasing speed, 
and simultaneously with the fracturing of the rock mass, the piezoelectric 
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voltage vanishes with the release of stress and the telluric currents instantly 
disappear. 

The telluric currents therefore occur and change as a precursor of an 
earthquake, and observation of the direction and intensity of the telluric 
5 currents therefore enables prediction of an earthquake. 

Since telluric currents do not flow on the earth's surface, direct 
detection of these currents is problematic, but induction fields that are 
caused by telluric currents (telluric current induction fields) do occur on the 
surface, and detection of the direction and intensity of magnetic fields on the 
10 surface therefore enables prediction of the direction and intensity of a telluric 
current induction field. 

The simplest method of observing a telluric current induction field is to 
observe the direction indicated by a magnetic needle. Under the influence of 
a telluric current induction field, a magnetic needle indicates a direction that 
15 diverges from the direction of normal terrestrial magnetism, and in an 

environment that is free of magnetism other than terrestrial magnetism and 
the telluric current induction field, a magnetic needle therefore enables the 
simplest and most economical observation. 

A second method allows more precise observation than a magnetic 
20 needle through the use of a magnetic force line sensor. 

A third method enables even greater accuracy by combining the use 
of a magnetic force line sensor and a GPS position detector. 

Alternatively, the combined use of a magnetic needle and a GPS 
position detector can obtain effective data for earthquake prediction. North 
25 as indicated by a magnetic needle does not align with true north, and it is 
further well known that this magnetic north shifts slightly every year. Thus, 
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true north is found by means of a GPS satellite, and through constant 
observation of the difference from north as indicated by a magnetic needle, 
effective data for earthquake prediction can be obtained. 

However, in order to raise the accuracy of observation data of telluric 
5 current induction fields, the magnetic field noise component that is 

generated by factors other than telluric current must be eliminated from the 
observation data. 

Magnetic field noise that occurs due to factors other than telluric 
current principally arises from the following sources: 

10 Magnetic fields that are generated by the direct current that flows in 

overhead power lines at observation points that are close to railway lines. 
Changes in such magnetic fields are characteristically short-term micro- 
changes, with the intensity of the magnetic field increasing as a train 
approaches and decreasing as a train moves away. 

15 Geomagnetic storms caused by the Dellinger phenomenon that 

accompanies solar activity. These disturbances characteristically occur and 
disappear in short time intervals. 

Magnetic fields that occur due to subterranean metal ores. These 
fields are characterized by constant fixed levels. 

20 In order to eliminate the magnetic field noise component that is 

produced by such factors other than telluric current, the magnetic field can 
be observed at fixed intervals at fixed observation points, the characteristics 
of the change patterns in these magnetic fields analyzed, and the magnetic 
field noise component then extracted and removed by software. 

25 FIG. 3 shows the earthquake prediction operations of earthquake 

prediction center 4. 
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First, telluric current induction field estimation unit 43: removes the 
magnetic field noise component at the observation point from observed 
magnetic field data (Step 101); then, as shown in FIG. 4, finds the amount of 
divergence from true north of the direction of the magnetic field at the 
5 observation point, from which the magnetic field noise component has been 
eliminated (Step 102); and estimates telluric current induction field N2 based 
on the vector difference between observed magnetic field N1 from which the 
magnetic field noise component has been eliminated and geomagnetic 
vector N that has been corrected to true north (Step 103). Then, as shown in 

10 FIG. 5, the results are plotted on a map (Step 104). 

Telluric current estimation unit 44 next connects the points on a map 
at which abnormalities in terrestrial magnetism have been recognized as 
shown in FIG. 5, and further, estimates the telluric current based on the 
right-handed screw rule (Step 105). 

15 Telluric current induction field intensity change pattern generation unit 

45 collects past data of the telluric current induction field intensities at 
specific observation points and generates a telluric current induction field 
intensity change pattern that shows change over time (Step 106). 

Earthquake prediction unit 46 analyzes the telluric current induction 

20 field intensity change patterns and the distribution of telluric currents that 
have been estimated by telluric current estimation unit 44, searches for 
unnatural regions such as points where the telluric currents are concentrated, 
and thus estimates seismofocal zones. Earthquake prediction unit 46 further 
compares and collates the telluric current induction field intensity change 

25 pattern that has been generated with past telluric current induction field 
intensity change patterns to estimate the time of occurrence and seismic 
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intensity of an earthquake (Step 107). 

As an example, when telluric current induction fields and telluric 
currents are plotted on a map of an observation area as shown in FIG. 6, the 
seismofocal zone can be estimated to be at point at which the telluric 
5 currents within the observation area are concentrated, and a large- 
magnitude earthquake can be estimated at a shallow layer directly below the 
observation area. 

Further, when the telluric current induction fields and telluric currents 
are plotted on a map of the observation area as shown in FIG. 7, the 
10 seismofocal zone can be estimated at a shallow layer in a vicinity outside the 
observation area. 

When the telluric current induction fields and telluric currents are 
plotted on a map of the observation area as shown in FIG. 8, the 
seismofocal zone can be estimated in a remote shallow layer outside the 
15 observation area, or in a plurality of close shallow layer points. 

FIG. 9 shows an example of the pattern of change in telluric current 
induction field intensity. In this case, the telluric current induction field 
intensity is a relative value. 

Typically, at a point close to the limit of elasticity that immediately 
20 precedes the breakdown of the rock mass of a seismofocal zone, the telluric 
current rises rapidly with the rapid rise in piezoelectric voltage. A stagnation 
of telluric current is then observed that accompanies the piezoelectric 
voltage saturation that immediately precedes plastic deformation, following 
which, the piezoelectric voltage vanishes with the release of pressure that is 
25 simultaneous with the breakdown of the rock mass, and the telluric current 
also disappears instantaneously. 



In addition, the transition in telluric current over time is specified only 
by the plasticity of the rock mass and relative vector speed of the two plates 
and has no relation to the distance from the observation points to the 
seismofocal zone; and since the telluric current induction fields are 
5 generated by telluric currents, observing the transition in telluric current 
induction field intensity over time by fixed-point observation enables an 
estimation of the transition of telluric currents over time. 

In this case, storing past telluric current induction field intensity 
change patterns of the environs of a seismofocal zone that occurs in a 

10 specific identical plate boundary plane enables extraction of telluric current 
induction field intensity change patterns up to immediately preceding the 
plastic deformation (occurrence of an earthquake) of the seismofocal zone. 
Accordingly, if the plates of an estimated seismofocal zone that is under 
observation can be specified, comparison and collation with past telluric 

15 current induction field intensity change patterns enables the estimation of an 
estimated time and seismic intensity until the plastic deformation of the rock 
mass (occurrence of the earthquake). 

In addition, transition points of a curve function that indicates change 
in the telluric current induction field intensity can be set and the time until the 

20 plastic deformation of the rock mass (earthquake occurrence) then estimated 
based on the telluric current induction field intensity at a point close to the 
limit of elasticity of the rock mass. Further, the maximum value that is 
reached by the telluric current induction field intensity can be estimated and 
the equivalent seismic intensity then estimated according to this maximum 

25 value. 

In the foregoing explanation, magnetic force line sensors are installed 
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in vehicles, ships, and existing fixed structures and magnetic field data are 
collected for each point within the observation area, but as another working 
example, magnetic force line sensors and GPS position detectors may be 
incorporated into mobile telephones or portable terminals, which may then 
5 use their own communication capabilities to transmit observation data. In 
such a case, the observation data may be periodically transmitted 
automatically, and if transmission fees are not charged, observation data 
can be collected from many points over a wide range without imposing any 
financial burden on users. 

10 Alternatively, if acceleration sensors are installed together with 

magnetic force line sensors on existing fixed structures and observation data 
then transmitted automatically when the acceleration sensors detect 
earthquake motion, magnetic field data can be collected at the time of 
occurrence of foreshocks that precede a main shock, whereby data effective 

15 for earthquake prediction can be obtained. 

Still further, if acceleration sensors together with magnetic force line 
sensors are incorporated in mobile units such as vehicles or ships as well as 
in mobile telephones or portable terminals, observation data may be 
automatically transmitted when the acceleration sensors detect a stationary 

20 state of at least a fixed time interval. 

As described in the foregoing explanation, earthquakes can be 
accurately predicted by collecting magnetic field data at a multiplicity of 
points within an observation area, estimating telluric current induction fields, 
estimating telluric currents based on the telluric current induction fields that 

25 have been estimated, and then analyzing these estimation results. 
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